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We present the analytic forms for the spectra of the cosmological perturbatfons 
from an initially anisotropic universe for the high momentum modes in the context 



of WKB approximations, as the continuation of the work [29[. We consider the Ein- 
stein gravity coupled to a hght scalar field. We then assume that the scalar field 
has the zero velocity initially and then slowly rolls down on the potential toward 
the origin. In the slow-roll approximations, the Kasner-de Sitter universe with a 
planar symmetry is a good approximation as the background evolution. Quanti- 
zation of the perturbations in the adiabatic vacuum, which we call the anisotropic 
vacuum, is carried out. For non-planar high momentum modes whose comoving 
momentum component orthogonal to the plane is bigger than the Hubble parameter 
at the inflationary phase, the WKB approximation is valid for the whole stage of 
the isotropization. On the other hand, the planar modes whose comoving momen- 
tum component orthogonal to the plane is comparable to the Hubble parameter, is 
amplified during the process of the anisotropic expansion. In the final gravitational 
wave spectra, we find that there is an asymmetry between the two polarizations of 
the gravitational wave because the initial mode mixing does not vanish. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Recent measurements by the WMAP satellite have suggested that the observed map 
of cosmic microwave background (CMB) anisotropy is almost consistent with the Gaussian 
and statistically isotropic primordial fluctuations from inflation, issues on a few anomalies 
in the CMB temperature map on large angular scales found in the recent data have been 
controversial. The most well-known fact is that there seems to be the suppression of the 
observed power of CMB anisotropy on angular scales bigger than sixty degrees [3| . There are 
other observational facts that imply the effect which induces the violation of the rotational 
invariance. More precisely, the planarity of lower multipole moments, the alignment between 
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the quadrupole (£ = 2) and the octopole = 3), and the ahgnment of them with the equinox 
and the echptic plane 0] were announced. There are other observational facts implying the 
large-scale anisotropy, i.e., odd correlations of £ = 4 ~ 8 multipoles with £ = 2,3 multipoles 
[1], a very large, possibly non-Gaussian cold spot in 10 degree scale j6[, asymmetry of angular 
map measured in north and south hemispheres [3] (see also the recent review [11 and more 
recent references therein). It also should be noted that some authors claim that there is no 
significant evidence for primordial isotropy breaking in five-year WMAP data joj (see also 
more recent arguments llO|). Indeed, to explain the origin of the anomalies, various solutions 
have been suggested, introducing a nontrivial topology (Tlj , a local anisotropy based on the 



Bianchi type Vllh universe to explain the quadrupole/octopole planarity and alignment [12 



non-linear inhomogeneities [13| and assuming an elliptic universe to explain the suppression 
of the quadrupole CMB power [3]. More recently, in particular, models which introduce 
an explicit source to break the spatial isotropy, either during inflation or in the late time 
universe, have been proposed, e.g., by the dynamics of an anisotropic energy-momentum 
component during inflation |l5l-ll8|. by the large scale magnetic field 19|, by the anisotropic 



cosmological constant [20( or dark energy [21[. 

The first purpose of this paper is to proceed to investigate the possibility that such large 
scale anomalies are produced by preinfiationary anisotropy and obtain the leading order 
corrections to the spectra for them. Cosmic nohair theorem ensures that in the presence of 
a positive cosmological constant an initially anisotropic universe exponentially approaches 
the de Sitter spacetime at the later time under the strong or dominant energy condition j22| . 
Therefore, it is plausible that the initial universe is highly anisotropic. 

The future CMB measurements will detect the fluctuations of B mode polarization in 
CMB, which may contain the information on the primordial gravitational waves. They 
would give a new tool to contrain the anisotropic cosmological model. The cosmological 



perturbation theory in the Kasner phase was formulated in Ref. 23l-l25l]. In general in an 
expanding (planar) Kasner phase one of two polarizations of gravitational waves is coupled 
with the scalar mode, but the other gravitational mode is decoupled. Thus, this coupling 
induces the asymmetry between propagations of two polarizations of the gravitational waves. 
If there are effects of the chiral symmetry breaking in the cosmic history, they would give rise 
to nonzero cross correlations between the fluctuations of the temperature and B mode, and 



of E and B modes [26|. They will give us powerful and independent tests on the primordial 
parity violation. The gravitational waves from the universe with an isotropy breaking would 
provide distinguishable signatures in the future CMB experiments. We will estimate how 
the initial mode mixing gives rise to the asymmetry between the primordial power spectra of 
the two gravitational wave modes, although we will not go into details of the observational 
aspects. Note that such subjects have been argued in the context of the anisotropic inflation 



models in Ref. [27|. The investigation of the higher order correlations as the bispectrum 



would provide us another interesting prediction to examine the the anisotropic universe (see 



e.g., [28|). We briefly comment on some expectations on this point in the last section. 

In the isotropic case, the quantization of fluctuations is carried out well inside the Hubble 
horizon, where the effects of the cosmic expansion can be ignored. In order to compare with 
the standard prediction, it is natural to quantize field in the initial adiabatic vacuum, which 
we call the anisotropic vacuum. There are two branches of the expanding Kasner solution 
with the planar symmetry. The initial adiabatic vacuum present only on one of the two 
branches where the expansion rate along the planar directions vanishes while that along the 
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special axis is finite ^. As a result, the initial spacetime structure can be seen as the prod- 
uct of two-dimensional Milne spacetime and two-dimensional Euclidean space. The scalar 
fluctuations decouple from the tensor fluctuations at the very initial time, so that the initial 
dynamics reduces to that in a system composed of three independent harmonic oscillators. 
In the other branch there is an initial singularity. Since the coupling diverges at the initial 
time, we cannot find an adiabatic vacuum. Therefore, here we focus on the first branch. For 
a given set of initial conditions, the power spectrum was investigated, rather by the numeri- 
cal ways in Ref. The aim of our study is to obtain more analytic understandings on 



the spectra from an initially anisotropic universe. Our previous work discussed the spectrum 
of a massless scalar field, ignoring its coupling with the metric perturbations [29[. In this 
work, as the continuation, we will discuss the metric perturbations, in particular focusing 
on the importance of the tensor-scalar coupling. 

The paper is constructed as follows: In Sec. II, the background solution of our anisotropic 
model is introduced. In Sec. Ill, we present the formulation of the coupled perturbations in 
the background of Kasner de Sitter solution with the planar symmetries and their relation 
to the cosmic observables. In Sec. IV, we investigate the behaviors of the perturbations 
modes after setting initial conditions in the anisotropic vacuum. In Sec. V, we close the 
article after giving a brief summary. 

II. BACKGROUND 

We consider the Einstein gravity minimally coupled to a massive scalar field 

S = l d'xV^[^R - ^-id<pr - \m^^) , (1) 

where g^j_i, is the spacetime metric and is a canonical scalar field with mass m. We consider 
an anisotropic spacetime with a two-dimensional planar symmetry 

ds^ = -dr^ + airfdx^ + h{Tf{dy^ + dz^), (2) 

where a(r) and 6(r) are independent scale factors. We define the expansion rates by 

Ha ■= -, Hb := -, (3) 
a 

respectively. In this paper, "dot" denotes the derivative with respect to the proper time r. 
Then, the field equations are given by 

H + 3H^='^, 3H'-h' = ^(l:^' + ^), + 3if0 + m^0 = O, (4) 



2M2 ' M2V2 



where the total and relative expansion rates are defined as 



Ha + 2Hb Ha — Hb 



^ The anisotropic vacuum is not specific to the Bianchi I model. In Ref. [29[, it has been shown that an 
anisotropic vacuum can also be defined for a Bianchi IX model. 
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Note that the evolution of the scalar field does not depend on the anisotropic scales but 
depends on the averaged scale factor. 

We assume that the scalar field stays initially at = 0o > and then slowly rolls 
down. Under this approximation, we can ignore the kinetic energy of the scalar field in the 
second equation of (jl]), and ~ 0o- From the first and then second equations of (HI), we 
obtain 

i/ = i/ocoth(3//or), h = ± . (6) 

smh(3iioT) 

where 

Ho = (7) 

is the total Hubble parameter at the asymptotic region in the limit of r — t- oo. Here, (-)- 
branch leads to the solution which contains a curvature singularity at the initial time which 
will not be dealt with in this paper. For the (+)-branch, we obtain 

Ha = . ,5° s h + cosh(3i7or)) , H, = Ho tanh (^H^t) . (8) 

At the initial time, they reduce to Ha — )■ \ and Hi, — )■ 0, which represents a (Milne) patch of 
the Minkowski spacetime and does not contain an initial singularity. The two-independent 
scale factors are given by 

a = smh^/\3HoT) tanhi (^^o^) , b = smh^^^SHor) coth^ (^^o^) • (9) 
The averaged scale factor is given by 

e" := (ab^) i = ( sinh(3iJo^)) ^ • (10) 

On this background, the evolution of the scalar field is approximated well by the behavior 
of a massive field in the above background 



v/?tan(f)Vr(l-|) 



r(l + |) ' '^2' 2 '^2' V' 



(11) 



where we defined x := GHqt, q := — IGM^ and M := = For q to be real, 

00 > \l\^p ""^ < 1^ the anisotropic stage Eq. (fTTj) can be further approximated 
as 

0^0o(l-^mV). (12) 

At the late time after the isotropization, r ^ l//Jo (2; ^ 1), oc g~2'f^o'^(i~^/i~i^^) . As 
long as 

00 > W -Mp, (13) 
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namely, m ^ the slow- roll condition is satisfied. 

Let us check the consistency of our approximation. The typical time scale for the cosmic 
isotropization is given by Xiso = 1, namely riso = On the other hand, the kinetic energy 
of the scalar field can be comparable to the potential at the time scale = ^. Therefore, 
the condition that </> almost stays around 0o and the kinetic energy of the scalar field is 
negligible is given by ^ riso, namely 00 ^ ^ from Eq. ([71). Thus, it turns out that with 
the slow- roll condition during inflation Eq. ( 1T3|) . under the assumption of the super-Planck 
initial amplitude, the solution Eq. (jH]) becomes a good approximation for the background 
evolution of the preinfiationary universe. 



III. PERTURBATIONS 



A. Perturbations in the Bianchi-type I universe 



1. Mode decompositions 



We then consider the cosmological perturbation theory in the anisotropic spacetime with 



a planar symmetry, following Ref. [2J]. The decomposition into the independent modes is 
performed in terms of the two-dimensional (y, 2;)-plane. A scalar quantity contains 1 degree 
of freedom. A vector quantity which satisfies the transverse condition in two dimensions 
also contains 1 degree of freedom. But there is no degree of freedom for a tensor quantity 
satisfying the transverse-traceless condition in two dimensions. Under the above decompo- 
sition, the totally 10 components of the metric perturbations can be classified into 7 scalar 
and 3 vector modes. Then, the perturbed metric can be described by 



( 



\ 



2$) 



adxx 
(1 - 2^) 



ad^B 



62fl-2S 



b'd,B, 
b^d2E3 



(14) 



where the matric is symmetric. £"3, S3 and B^ correspond to vector modes, and the rest 
are scalar ones. In addition, there is the perturbation of the scalar field, -f 50, which is 
definitively the scalar mode. Thus, there are totally 8 scalar and 3 vector modes. 

In the rest of the paper, we work in the momentum space after decomposing perturbations 
into (comoving) Fourier modes. We distinguish the comoving momenta ki from physical 
momenta pj, which depends on time, as 



Pi 



a 



The total momenta are defined by 



P2 



b ■ 



(15) 



hj — k-^ ~\~ k' 



25 



p 



2 I 2 
Pl+P2- 



(16) 



Thanks to the residual planar symmetry on the (y, z) plane, we can fix A;3 = without any 
loss of generality, ki denotes the component of the momentum along the special x direction, 
and k2 does that in the orthogonal plane. 
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2. Vector mode 

About the vector mode, 1 of 3 components can be eliminated by the gauge fixing and the 
other 1 can be done by the momentum constraint. By setting the gauge of vector mode to 
be E's = and ehminating the nondynamical component -B3, we have only one propagating 
normalized degree of freedom 



V2 



i,2 I a^^2 



-.Br, 



(17) 



whose equation of motion is given by 



+ »^ 



with 



pl+pl-Ha{Ha-H,)+H, 



2M2 



[Ha - H,) 



(18) 



(19) 



For convenience, we introduced a new time coordinate t by 

dr dr 
dt = — = — 



ah' 



3a 



(20) 



in which 



sinh ( 3ifoT 



sinh (3Ho{-t) 



Note that t —00 as r — )■ 0, and t — ?■ 0— as r — )■ 00. 



3. Scalar mode 



About the scalar mode, 3 of totally 8 components can be eliminated by the gauge fixing 
and the other 3 can be done by the constraints (1 Hamiltonian and 2 momentum constraints) 
and hence there are two propagating degrees of freedom. Setting B = Ti = E = Ohj fixing 
the gauge, and then eliminatng the nondynamical components $, x ^'^'^ B, two normalized 
propagating degrees of freedom are given by 



V 



y/2MpplHb 



" H,pI + H,{2pI+pI) 
V and obey the coupled equations of motion 



H,pI + H,{2pI+pI) 



■ d' 



^11 ^12 



^12 ^22 



V 



(21) 



(22) 
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where the components of the frequency matrix are given by 



11 



2M2 



+ 



2H„ 



22 



a'hH pi+pi + Ht{Hh - Ha) + Ht + 



2M2 



pl{Ha-H,y 



AHu 



2H,pI + {Ha + H,)pI 



12 



3^+ 1 



2Hy, + (Ha + H,)pI 
, , V2pl{Ha - H,) 
2H,pI + (Ha + H,)pI 

pliHg - H,) / 1 



+ 



2 



/f, M2 2HiM^ Hb Ml 



+ m 



2H,pI + + H,)pI Mp L H, H^M^ H^Mp 

pljHa - H,) ' 
2H,pI + (/f, + Mp 



2i7. M3 M„ 



(23) 



4- The isotropic limit 

Before closing this section, we briefly mention the limit to the ordinary homogeneous and 
isotropic universe 6 — )■ a. In the later times, we find 

\ Mm 2M2 2MmV ^ z ° ^' 



22) 



^ «^ f^' + ^) = «^ f - - + '^Hy) , (24) 



2M2, 



a 



where z = while — )■ 0. The prime denotes the derivative with respect to the 
conformal time dr] = The late time evolutions are given in terms of those in the de 
Sitter spacetime, written in terms of the Bessel functions. Thus, in this limit, V reduces 
to the Sasaki- Mukhanov variable v of scalar perturbations in three dimensions, defined in 



30[, while H^ and H^ give two independent tensor polarizations and /i+ in the flat 
three-dimensional space. Therefore, through the propagations in the anisotropic universe, 
an asymmetry between the two tensor polarizations would appear. 



B. Quantization and power spectra 

After giving the equations of motion, we are going to quantize the perturbation modes. 
At the initial times, the frequency squares of the vector and scalar modes in Eqs. (fT9|) and 
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( !23|) . behave as 



= 2I14 + ( 2HAkl + 3kl) + ^ + 0{r\ 



kl 



.3n 



108HM \ 3Hh^ 



ul, = 2tfc? + ( 2H^kl + 3kl) + ^f^^ ] ^ + 0(r^^ 



2 36v^ifom 3 

= - 4fc? + 3kl ^ + ^^^)- (25) 



The frequency squared appears to be divergent in the hmit ki — 0. However, it is 
an artifact of the early time hmit as one can see in Eq. (fT9|) . Since all of uj\, loI^ and 

approach constants, and — as r — )■ 0, namely, the coupling between V and 
vanishes in the early time, the adiabatic vacuum can be found. To distinguish this vacuum 
state from the standard Bunch-Davis vacuum, we call our adiabatic vacuum an anisotropic 
vacuum. Then, these perturbations are quantized in this vacuum. The procedure of the 
canonical quantization follows the standard manner. 

Here, we introduce Y for the collective notation of i^x? V and For the anisotropic 
vacuum |0), the anihilation operator can be defined by ak|0) = 0. Then, the metric pertur- 
bations are canonically quantized as 

Y = J d^k(ukai, + ulal^, (26) 

where operators satisfy the commuation relation [aki, a|^^] = (5(ki — k2) (others are zero) and 
= e'*'^Yk/(2vr)^/^ are the mode functions satisfing the normalization condition Yi^dtY^ 



{dtYk)Y* = i/e^". Frow now on, we omit the subscript k. Since from Eq. (1241) the mode 
mixing is also absent in the late time universe, the final power spectra of these modes are 
independently defined by 



P - Iff |2 



(27) 



and related to the late time behaviors of perturbations. For comparison, in the standard 
slow-roll inflation, the spectra are given by 



p(o) _ { p{0) _ p(o) _ 

~\27r^) ' ~ -2nm^- ^^^^ 



C. Evolutions of perturbations in the intermediate times and WKB approximation 

In the intermediate times, two modes are generically coupled. If the WKB approximation 
is valid, however, we may simplify to solve these perturbations. The original equation 
Eq. (!22l) can be described by the matrix differential equation 



(29) 
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where 



Defining M = AAi, where A is a rotational matrix 



A 

Eq. (l29i) becomes 



cos 6 — sin 6 
sin 6' cos ^ 



+ -JT-JT + U -37^ + n')M = 0. (31) 



df^ dt dt V dt"^ 

By choosing 6 so that Cl"^ = A"^Q'^A becomes a diagonal matrix, we obtain 



tan(2e) = -X, tan^= -±./^ + l, X := 2^"^'" 2 • (^2) 

X V X "^22 ~ ^11 



COS 6* 



1 + 



sin 6 = ± 



1 - 



(33) 



•^22 > '^ll 



^^22 < '^ll; 



where the positive (negative) sign will be taken for cu. 
the limit x — ^ =F0. Note that x < corresponds to 6* > and vise versa. 
At the initial period of time, 



X 



2m/ Ho 



V3 



1 + 



6H; 



+ 0(t' 



to have 6' = in 



(34) 



is of 0{m/HQ). Thus, interestingly, although uji2 — as t — ?■ — oo, x remains finite in the 
same limit and there is a finite amount of the mode mixing initially. On the other hand, for 
t — )■ oo, X 0. The effect of the mode mixing is encoded into terms which contain the time 
derivatives of the rotational matrix A in Eq. (I3T!) . We also find 



A-^n^A 



~,2 





ul^ 



11 



where 



,~.2 




(^^22 



UJ 



22 



^11 + ^^22 + 




(35) 



Then, 



M 



V 



V" cos 6* + if + sin 6 
if _i_ cos 6' — y sin 6' 



(36) 
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where ^ is a slowly varying function as 9' 
valid. Thus, 



0{uj[i^uj'22) SO that the WKB approximation is 



\uj^M\ > 



A 



-1 



A]—M 



A 



dt J dt 



dt 



ujM 



-M 



Therefore, to obtain the leading order correction, it is enough to solve 

d'^M 



df^ 



+ n'M ^ 0. 



(37) 



(38) 



If the WKB is not valid, one has to solve Eq. (j^HD directly. 

Introducing a new corrective notation Y for , V and then we can check the validity 
of the WKB approximation through the adiabaticity parameter. 



where represents either of 



ey 



I i2 r,2 



1141 



(39) 



or Cj22- Foi' the high momentum case, since the 



behaviors of Uy are very similar to to the case of a massless scalar field discussed in \2{ 
from Eq. ( 12^ . e-p also behaves as in the same way (and see Fig. I). Hence, for the non-planar 
high- momentum modes ki ^ k2 ^ Hq, ey is always less than unity in the early times, while 
for the planar modes Hq < ki <^ k2 it temporarily exceeds unity, which implies that the 
WKB approximation is broken there. Therefore, in the next sections we treat these two 
cases separately. 



IV. NON-PLANAR MODES WITH A;i, /cz > ^^o 

For the non-planar high momentum mode, the WKB approximation is valid for a small 
scale factor. We can discuss the evolution of the perturbation modes by solving Eq. (1551) . 
Since the equations are approximately diagonalized, we solve them correctively. 

In the early time, all H-^ , V and if+ modes are quantized in the adiabatic vacuum. Their 
initial amplitudes are given by 

^xL =y\, =H+\, =-^. (40) 

On the scalar mode from Eq. f l5^ . at the initial time the mixing angle is not vanishing, 
sin^ tIh" ^^'^ cos^ ~ 1, hence initial amplitudes of modes after the diagonalization are 
given by 

V\ ^{v ^-!^e\\ , ^(E^--!^y\\ , (41) 

i— !> — no \ /orr ' / — no' ^ — oo \ ^ fn tt I !> — no' ^ ' 



respectively. For the high momentum modes, ignoring the corrections of O(m^), = iiJ22 = 
a)ii = bJ22 which are correctively represented by u. In addition, since the WKB approxima- 
tion is valid for the non-planar high momentum modes, we obtain the time evolution in the 
intermediate times as 

V^^(l + ^^)^e-^/*'^*'^, iJ+^(l--^)^e-^^^'^*'^, (42) 
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where Q is given by 



2\ n 



3^ 
2 1]2 



and satisfies Q ^ ki for t — )■ — oo. On tlie otlier liand, 

1 r* M'( 



i ydt'n 



(43) 



(44) 



By using b ~ (— 3iJo^) ^ and —HqI] ^ {—3 Hot) I , since tlie mixing is eventually vanislies, 
tfie late time solutions are given in terms of the de Sitter mode solutions by 



Y 



2k 



B 



Y 



;45) 



We now match the WKB solution to the de Sitter solutions. The details of the matching 
are summarized in Appendix A. Keeping the accurary to the second order in adiabatic 
approximation, we finally obtain the power spectrum including the leading order corrections 
from the direction dependent part 



Px 

Pv 

P^ 



P. 



(0) 



Hn\ 2 



k 



l + COS 2W-- +0 




k 



p 



(0) 
V 



p 



(0) 



1 + 
1 - 



2m 
2m 



Hn\2 



Q(r2) 



k 
k 



+ g(r,)(^) cos(2^/-)+0((^^ 




\ 2 



COS ( 2 



k 

Ho 




O 



k 



(46) 



2 -r| (ra : = 



— ) denotes the leading order corrections due to the anisotropy. 



where Q{r2) .-3 , ^ ^ z ^ 
Thus, the leading order correction appears at the third adiabatic order and mode mixing 
gives the contribution of order j^, which results in a chiral asymmetry between two tensor 
polarizations. 



V. PLANAR MODES Hq < hi ^ fca 

For the planar modes of ki <^ k2, the WKB approximation is violated in the early times. 
Instead, we make use of approximate solutions. We expand the effective frequency as. 



where 



2 



2I {kl + e^^"*) , oo^, ^ 2I {kj + Soo^.e''''') , 
'.2l{kl + 5ul,e''''^'') , Lof, ^ 2I (5a;2^)e6^»*. (47) 



4 Qff2l,2 4 41.2 _ ol2 

_ 12 . ^,2 I ^-"0^2 r 2 _ 12 ^,2 o2/3 2 ^«^1 

'^^x " ^2 + 3^^1+ 2i^2 ' ^^11-^2 + 3^:1 2 ^^4^2 + 3^2' 

4 18x2^3^ 6 X 2^I^V^Homkl 

5u^,, -k, + -k, + 4^2 + 3^2 ' ^^12 - ^^^2 • (48) 
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For convenience, we introduce the planarity parameter s 1), defined by /c2 = 
Then, 6us can be expanded around s = as 

23 , 2 , ^/ 2n r ,2 _ '^1 , 



Slo 



kl /4 



'^^22 = ^ + (3^1 + 6 X 2§i7o') + 0(5^), <5w?2 = -2lV3Hom + 0(5^), (49) 



and 



' ' A;? + 23 x3ifo%l-^/l + ;^ + ^)+2-3m2 + 0(.2). 



^^2 _ kl , 4,2 , ,1 ^^2^ , /i , , 



" s2 ^ 3 



fc? + 2§ X 3i/o^(l + Jl + ^ + ^) + + 0{s^). (50) 



'0 ^W-I-Ig. 

In the hmit of s — )■ the frequency diverges and we focus on s 7^ 0. Note that 



2 / TTi^ Tfl^ 2 

6Ql, - 6Ql^ = -6 X 2^HlM+ + ^ ^ -6 x 2^ Hi < 0. 
In the small s limit, a;|, behave as 

4 = 2ifc^inht ( ^) sinhi (3i7or) + 0(s°) = + 0(s°), (51) 

where the species of perturbations are irrelevant. 



A. Hx mode and power spectrum 

For Hot ^ — 1, the solutions for Hx are given by 



= X„ " (gxe^^°*) , (52) 

2 23(50;^ 

where gi := and qx '■= — 3^^^ — • Note that this solution reproduces the correct normal- 
ization in the limit t — )■ — 00. 

When the WKB solution becomes valid, we can match the solution Eq. (!52l) to the WKB 
solution. Then, when the universe enters into the de Sitter phase, we can match the WKB 
solution to the de Sitter mode functions. The details are summarized in Appendix B 1. In 
this subsection, we show the final power spectrum of Hx for the planar modes 

where 

v?r(i) k rr 

Here we use the fact that ^2 — /c in the case of planar modes. 

As shown in Fig. [H the power spectrum exponentially approach to that in the isotropic 
case. On the other hand, for the planar mode with ki/H^ be 0(1), it leaves observable effect 
in the CMB. 
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FIG. 1: Power spectrum of the tensor mode Hx relative to the isotropic case. The high momentum 
mode corresponds to large ki/HQ and k2/Ho. The planar mode corresponds to the case that ki/HQ 
becomes 0(1). 



B. Mixed modes and power spectra 



The evolution of V and if+ modes can be dealt with the parallel way. The evolution 
equations at the initial times are given in terms of the mixed form 



V 



24/3 



kll 



V 



0, 



(54) 



where 5u'^s are defined in Eq. (148|) . We may diagonalize the equation (154|) by using the 
transformation 



By choosing ip to be 

COS'?/' = 



V 



rl 



O 



cos ijj — sin ip 
sin ip cos ip 



xl 



At the initial times, xo becomes 

Xo = 



2m/ Hq 



v^(l-mV(6i72)) 
Then the frequency term can be diagonalized 



+ 0(1/5^ 



o 



-1 



O 



5Coi^ 
5Cj^ 



22 



(55) 
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where 



2 V 



{5, 



22 




^^22 




A{5UJ 



12) 



11 



^ {dul^ + 5^22 + (^^22 - ^^^ll) 

The original mixed equations are now diagonahzed for the new variables 

M = 



(56) 



V cos ip + sin ip 
cos ip — V sin ip 



(57) 



The solutions are given by 



m 



m 



TV 



QHq sinh(7rgi 



-X 



TT 



VSHoyy QHq sinh(7rgi 



-J_ 



iiji 



que 



3Hot 



q22e 



3Hot\ 



(58) 



where qi :- 



2i|fci| 



- ---2 ^2 



23 (5(1, J 
3^ 



2TJ(5i>2^)7 



and similarly to the non-planar 



3Ho ' W ' ^22 •- 

high momentum modes, factors of order represent corrections due to the initial mixing of 
modes and the early time limit of the above solutions correctly reproduces the normalized 
amplitude. Note that these solutions reproduce the correct normalizations in the limit 
t — )■ — oo. By definition, V and iJ+ are now given by 



V 



TT 



1 



m 



QHq sinh(7rgi 

m 



cos ipJ- 



que 



3Hot 



3Hot 



TT 



6i/osinh(7rgi) LV^ 



smipJ_igAque 



^226 



3Hot 



(59) 



When the WKB becomes valid, the coupling between V and is negligible. At this 
time, we can match the solution Eq. flM?]) to WKB solutions. Then, when the universe 
enters into the de Sitter phase, we finally match the WKB solutions, to the de Sitter mode 
functions. The essential procedure is almost the same as the case of the mode. The 
details are summarized in the Appendix B 2. In this subsection we only show the final 
results: The power spectrum of V (= V^) for the planar modes is given by 



V 



(0) 



2m 



where 



VA 



v^r(| 



sinh(7rgi) / ' 



k 



(60) 



3 X 2^r(§)^o 
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Similarly, the power spectrum of iJ+ (= for the planar modes is given by 

- V VSHo) y ^ sinh(vrgi) ) ' 



3 X 23r(|) -no V -no 

Note that the forms of the power spectra for the and V modes are almost the same as 
that of the mode except for the global scaling due to the initial mode mixing. 



VI. CONCLUSION 

In this paper, we have analytically investigated the corrections to the power spectra of the 
cosmological perturbations due to the preinfiationary anisotropy of the universe, for the high 
momentum modes in the context of WKB approximations. The first motivation to consider 
the anisotropic universe is that even if the present universe is almost isotropic, it does not 
mean that it is also isotropic from the beginning. It would be more generic that the initial 
universe is highly anisotropic. The second motivation comes from observations. In recent 
years, several groups have reported the so-called low-£ anomalies in large angular power of 
CMB fluctuations. They may be produced by the breaking of the rotational symmetry in 
the early universe. 

We considered the Einstein gravity coupled to a light scalar field. We assumed that 
this scalar field initially stays at a very large (super-Planck) field value, and then starts to 
roll down slowly. If the mass of the scalar field is small enough, the kinetic energy of the 
scalar field does not affect the spacetime dynamics significantly. Imposing the regularity 
of the spacetime at the initial time, one of two planar branches of the Kasner-de Sitter 
solution, whose initial geometry becomes a (Milne) patch of the Minkowski spacetime, is 
a good approximation for describing the cosmic isotropization. The cosmic isotropization 
takes place within a few Hubble times. During the subsequent inflationary stage, the scalar 
field rolls down toward the true minimum and plays the role of the inflaton. 

Then, we investigated the analytic expressions for the spectra of the cosmological per- 
turbations produced in the above anisotropic background for the high momentum modes. 
In the anisotropic background, there are 2 scalar and 1 vector modes in terms of the two- 
dimensional flat space. In the isotropic limit, this vector mode in the anisotropic universe 
reduces to one of tensor polarizations in the flat 3-dimensional space, while two scalar modes 
reduce to one scalar and the other tensor polarization there. During the anisotropic phase, 
two scalar modes are coupled, resulting in the asymmetry between the spectra of two ten- 
sor polarizations obtained in the isotropic limit. Since at the initial times, the coupling is 
absent, we could define the adiabatic vacuum, which we call the anisotropic vacuum, and 
canonically quantize the perturbations. Our anisotropic vacuum is definitely different from 
the standard Bunch-Davis vacuum. In addition, the presence of the anisotropic vacuum is 
specific to our particular choice of the Kasner parameter. 

For the non-planar modes k,ki ^ Hq, for the sufficiently high momentum the WKB 
approximation is valid but the mixing angle does not vanish at the early time. At the 
leading order, the power spectra of all the perturbation modes contain the corrections due 
to the nonstandard propagation in the anisotropic background, and the effects of the initial 
mixing of modes. The former has the universal form, while the latter induces an asymmetry 
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of two gravitational wave polarizations in the isotropic limit. The modifications of spectra 
appear in the oscillatory behaviors of the primordial spectrum. On the other hand, for the 
planar mode, i.e., k ^ ki ^ Hq, although the WKB approximation is broken at the very 
early times, the mode mixing does not take place significantly. For the modes of fci = 0, 
the adiabaticity parameter initially diverges and hence the anisotropic vacuum is not well 
defined. However, such modes are not relevant for the observations. 

One of the main results which we have obtained is that, irrespective of the non-planar or 
planar modes, for the high momentum modes the ratio of the power spectra between two 
tensor polarizations is given by 

^ ~ 1 - (62) 

In the chaotic inflation typically = O(O.l), and hence there is a difference of the spectra 
of two gravitational wave polarizations, which is of roughly ten percent. 

Before closing this article, it may be important to mention the effects of the primordial 
anisotropy on the higher order spectra, in particular on the bispectrum, and the primordial 
non-Gaussianities. In the model discussed in this paper the initial vacuum is not the standard 
Bunch-Davis vacuum, and the subsequent evolution is almost the same as that in the single 
field, slow-roll inflation. Thus, we expect that the folded shape bispectrum where three 



momenta satisfy ki + k2 ^ ks would become dominant |28l. l3l[| (see also |32|), while the local 
shape bispectrum of ^2 ~ ^3 2> fci would be negligible 33|. However, the particular local 
type bispectrum where and k^ almost lie in the plane of the y and z directions of Eq. ([2]) 
while ki is nearly orthogonal to this plane (^2 ~ ^3 3> fci), could be much different from the 
case of the isotropic universe, since as in the case of the spectra for the planar modes it could 
be sensitive to the anisotropy. For the non-planar case of ki + k2 ^ ks, the amount of the 
non-Gaussianities would be determined by the coefficient By in Eq. (H5l) which represents 
the amplitude of the negative frequency mode; 

'Hn\i . , rr 



ReiBy) ~ g(r2) (^) %in (^— + 0o) , (63) 



where 0o denotes some constant phase. Note that due to the effect of the initial mode mixing 
Re{By) may be amplified by some factor of order ~ 0.1. Therefore, the bispectrum could 
exhibit an oscillatory behavior, and contain the information on the primordial anisotropy 
through the factor Q{r2), which may be distinguishable if detected. The concrete evaluation 
of the bispectrum in the anisotropic universe and the detectability of the non-Gaussianities 
will be interesting issues and should be investigated in the future studies. 
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Appendix A: On the power spectra for the non-planar, high-momentum modes 

We will summarize the matching of the WKB and de Sitter mode functions for the non- 
planar high momentum modes, given by Eqs. ( l42l) and ( l45l) . respectively. 

Keeping the accuracy to the second order in adiabatic approximation, the solution to Eq. 
I3|) is given by 

n = uj{i + s), (Ai) 



where 



Defining 



4^2 V uj 2io^ 



and using the results of matching at some appropriate time t = t^, the final power spectrum 
is given by 

+ (-{a + c^- ^^)sin[^(-3i/ot)^] +ci^*cos[^(-3i/ot)^])'|,(A3) 

where z/^, := ^(— 3ifo^*)^ and $ = /** dt'ui{t'), and Py^ is given by (!28|) . i^y represents the 
contribution of the initial mixing of modes, which are given by = 1, Ky = 1 + and 

^ _ 1 2m 



For all and modes, we obtain in the high momentum limit, 

a;2 = pg4a _^ 2e-=^°Q(r2)) + 0(m2), (A4) 

where Q{r2) '■= \ —'^2 {^2 •= ^) denotes the leading order corrections due to the anisotropy. 
Note that — | < Q{r2) < | and Oiw?) terms contain the slow-roll corrections. Thus, as long 
as our concerns are in the high momentum modes, 0{rrP') terms are not important. The 
WKB approximations are valid as long as the adiabaticity parameter Eq ( 139|) is smaller than 

unity, namely —HqI ^ (^) • other hand, the late time approximation is valid if 

e" ^ 1. The matching time t = should be chosen so that the error is minimized and is 



given in terms of their geometric mean [29 



6°= 



-Y,* — 

is given by 



(A5) 



1 

At t = ey^^ - ^. Making use of {-3HoQi = (f^) ' 1 + 0(| 



■ 



. the power spectrum 



+ ( - (a + cz/.y^) sin + cos Y^)'}- (A6) 
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As shown in the previous section, sin^ is 0(771/ Hq) at early times and vanishes at late 
times, and thus the initial condition for V and modes will be specified by using the 
WKB approximation for V and Since V = V and = at t — )■ 0—, the power 
spectra for Y reproduce those of Y. Therefore, Py = Py- 

Then, evaluating at t = ti up to 



u, = l + g(r2) 



\ 2 



O 



k 



a 



2k \\ k 



2k \\ k 




Ho 

k 



Ho 
2k 



(A7) 



Therefore, we obtain the power spectrum including the leading order corrections from the 
direction dependent part 



Px 


_ p(0) 

X 


1 + 


Q{r2)( 


Pv 


_ p(0) 


1 + 


2m 
V^Ho 


P+ 


= Pi'^ 


1 - 


2m 
VSHo 



Hp \ 2 

T 




k 

Ho 



+ g(r,)(|^)%os(2,/- 




Q{r2) 



Ho 



cos ( 2 



k 

Ho 
k 




Hn 



O 

o 



Hoy 
k J 

Ho\^ 
k 



(A8) 



Appendix B: On the power spectra of the planar modes 

1. On Fx 

Here, we explain the derivation of the power spectrum of Hx for the planar modes. 
For Hot <^ — 1, the solutions for Hx are given by 



H. 



(1) 



IT 



6H0 sinh(7rgi 



-J. 



3Hot\ 



(Bl) 



where qi := and qx '■= Note that this solution reproduces the correct normal- 



ization in the limit t — t- —00. During the period of time 2ti x < t ^ where ti is the 

matching time given below, the WKB approximation 



H 



(2) 



B 



V2LJ, 



=exp 



dt'cux (f) 



l,x 



Bx^ 



=exp 



dt'uJx {t') 



(B2) 



holds well. In order for the error to be minimized, we need to choose the matching time of 
the two solutions (IBll) and (]B2p . t = ti x to be 



^0 

2§ V {5ulf2 



(B3) 
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since the adiabaticity parameter behaves as ex — ^^^^^° e ^^ot_ Since q^e^^°^^''' 1 at 
t = ti, 



Q rr + / \ Sin TC* ^ ^ 
SHoq^e-^^o* sinh(7rgi) 



(B4) 



where = gxC^^"* + f + Noting 



|t=ti 



2 X _l 

2 



(B5) 



we obtain 

-Bx -I- = 



rril) 



sinh(7rgi) 



5. 



l + -^r'cosT 

(1) 



X U=tl,> 



-2 1 



Hoqt 



2 ^ 1 



- sinT> 



it=ti, 



2 sinh(7rgi) 



1 H r I COS T 



(5u;^x) 



.2 , 1 



xU=t.x+^(l + -^)'sinTx|,^,^ 



Since <^ — at the leading order they reduce to 



-'-'x,+ — 



-Txqi 



(1-e- 



27rgi^ 



) -Dx,— — 



(1 _e-2^9i)^ 



(B7) 



Then, at t = ^ l^i.xl), where is chosen to be Eq. (lASp as for the non-planar 

modes, the WKB solution are matched to the de Sitter mode function 



H 



(3) 



'2k 



k 



For the high momentum modes, 



C 



2k 



( - SHot) 



k 
(B8) 



Cx,+ ^ e" 



"e-^*5x+, Cx- ^e'V^^oe'^S 



'X,+ ; >-^X 



X,-' 



(B9) 



where $ = //* dt'ux{t'). Thus, 



Cx,+ ~ Cx, 



(1 - e-2'f9i)2 



r<! (1 - e--"^) cos - ^x) + ^(1 + e-'^^^) sin - ^x) felO) 



vr 
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where \I/x 
by 



^^0 



Therefore, the power spectrum for the mode is given 



Ph. ^ <^(coth(7rgi) 
Finally, we estimate \1/ x • From Eq. (15T]) , 

12 3 



sin(2^x) \ 
sinh(7rgi) / 



dt = -r^e^^°\Fi 



;i_e6Hot)f 
and the $ term dominates \l/x: 



2' 3' 2'^ 



A: 




3 X 2ir(|)iyo^ 



(Bll) 



(B12) 



Here we use the fact that k2 — k in the case of planar modes. 



2. On V and H+ 

Similarly, we explain the derivation of the power spectra of V and iJ+ for the planar 
modes. Firstly, we match the solutions Eq. ( l59|) to the WKB solutions, and then match 
them to the mode functions in the de Sitter inflation. Most computations are the same as 
in the case of ■ 



a. V mode 

At the time when the WKB becomes valid, the coupling between V and is negligible. 
At t = tiv given by 



e3^^o*M'^ 4 /_^, (B13) 



the asymptotic form of the first solution 

\ cos [ip — 9^ 



m \ sin (ip — 6^) _ "I ^ 
v3-Ho^ 23(d)a;22)'* 

where Tn = gne^^o* + f + ^ and T22 = q22e^^°^ + f + is matched to the WKB 
solution given by 
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Since from Eqs. ( 133|) and (155|) . at t = ti 



X 



12 



2(5^2 



12 



Xo, 



(B16) 



at the leading order, we obtain ip ^ Oiy. The first solution then reduces to 



V 



Noting that q'f^ <ti ^'^'^i)^ ^ after matching these solutions, we obtain 



m 



— sin T 



11- 



(BIT) 



Ho 



Bi 



(1 - e-2'rgij2 



r 1 



(1 - 



91 ) 2 



(B18) 



Then, at t = ^ |^i,y|), where is chosen to be Eq. (lASp as for the non-planar 

high momentum modes, the WKB solution are matched to the de Sitter mode function 



'2k 



e "0 



+ 



2k 



e "0 



For the high momentum modes, 
1 



k 

(B19) 



Cv,+ — Cv,- 



^— ^{ (1 - e--0 cos - ^y,) + e--0 sin - ^v^l^O) 



where ^yi := \l + $ + and *y,2 : = 

spectrum for the V (= V) mode is given by 



- + $ + 1 / -^. Thus, the power 



Pv ^ pfd 



where the $ term dominates '^v,i and 

v^r(|) k 



2m \( sin(2^y,i)\ 

coth(7rgi) - — TT r 

sinh(7rgij / 



(B21) 




Here we use the fact that A;2 — A; in the case of planar modes. 
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b. H+ mode 



Similarly for the Hj^ mode, at t = ti^ 



2a y (dui^)^ 



the asymptotic form of the first solution 



1 



1 + ) 1 „ . 1 sm Til 



+ 1 



is matched to the WKB solution is given by 



sinh(7rgi) Lv \/?,Hq^ 2\{6ujl^)-i 

m xcosiip - 9) 

— cos T22 



H 



(2) 



\J2U}'22 



exp 



dt'u22{t') 



y/2u22 



exp 



dt'uJ22it' 



(B22) 



(B23) 



(B24) 



From Eqs. fl33|) and fl55l) . at t = ti^+ the similar relation to Eq. flB16P is obtained, which 
leads to ~ 6'i _|_. Thus, the first solution reduces to 



, 1 / m \ 1 „ 
~ 4 / , ■ . . rll ^ — I— ^ z — ^cosT22. 



23 {5Q. 



~2 ] 
22) 



(B25) 



Noting ^ ^~^o^' cifter matching these solutions, at the leading order we obtain the 
coefficients 



B, 



1 



[1 - e-2'r''ij2 



r 1 



m 



^^0 



) 



m 



(B26) 



Then, at t = with ^ where is chosen to be Eq. f lASp as for the non-planar 

high momentum modes, the WKB solution are matched to the de Sitter mode function 



H 



(3) 



2k 



( - SHot) 



*-^0" V2F V 
For the high momentum modes. 



k 

(B27) 



(1 - e--) cos - ^^^2) + sin - vl/_,,2|^?8) 
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where ^+ i := ,/ + $ + .ff and ^+ 2 := \/ ^^^^^ + $ + Thus, the power 

^' Y Ho(5tD22)i V ^0 V ^0 V ^0 

spectrum for the iJ+ (= mode is given by 

- y VSHoJy ^^'' sinh(vrgO J" 

v^r(|) k [T 



= . .il'' Tr + Q(\/Tr)- (B29) 



Here we use the fact that k2 — k in the case of planar modes. 
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